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ABSTRACT 


The primary electron spectrometer used to detect auroral electrons 
on sounding rocket 18:63 UE is described. The spectrometer used 
exponentially decaying positive and negative voltages applied to 
spherical deflection plates for energy analysis. A method for deter- 
mining the analyzer response which does not require the assumptions 
that the ratio of plate separation to mean radius, the entrance or 
the exit apertures are small is described. By comparison with 
experiment it is shown that the effect of neither entrance nor exit 
collimation can be ignored. The experimental and calculated values 
of the limiting orbits agree well. A non-iterative technique of un- 
folding the electron differential energy spectrum is described. This 
method does not require the usual assumption of a flat or histogram- 
type energy spectrum. The unfolded spectra using both this technique 
and one which assumes a flat spectrum are compared to actual input 
spectra. This technique is especially useful in analyzing peaked auroral 


electron energy spectra. 



I . Introduct ion 

t 

Theocloridis and PaolittL have reviewed the problem of determining 

2 

the differential particle flux (particles/sec-cm -sr-lceV) from the 

number of counts accumulated during a given time interval with fixed 

symmetric voltages on the deflection plates of an electrostatic analyzer. 

2 

In a previous paper they described a method for calculating an energy- 
angle factor for spherical plate electrostatic analyzers. This energy 
angle factor which they call <AotAE> describes the particle transmission 
in the plane of the particle trajectory as it passes between the analyzer 
plates. Because of the coupling between energy and angle it is the 
motion in this plane which makes the determination of the response to 
an isotropic flux for an electrostatic analyzer differ from that of an 
analyzer where the geometric factor is decoupled from the energy deter- 
mination as for example a collimated solid state detector. With the 
usual assumptions of a central force analyzing field and neglecting 
fringing fields the response as a function of the orientation of the 
plane of the trajectory is decoupled from the particle energy and can 

be determined from purely geometrical considerations.^ 

3 

Smith and Day have reported an improvement upon the Theodoridis 
and Paolini technique of determining the analyzer response which does not 
require the assumptions that the ratio of plate separation to mean radius, 
AR/R, is small, that one of the angular acceptance windows Aa (azimuthal) 
or A3 (polar) is small, and that the energy-a response is independent of 
the polar angle 3. The method of determining the analyzer response function 
which we use is an improvement upon the Smith and Day technique in that 
the effects of entrance and exit collimaticn are included in the determinate 
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of allowed orbits. In the electrostatic analyzer we describe neglecting 
the colliination effects can produce a factor of two error in deter- 
mining the differential flux. 

A. Detector Description 

An electrostatic analyzer with partial spherical geometry was used 
for electron energy discrimination in the Primary Electron Spectrometer 
(PESPEC) used to measure auroral electrons on Nike-Tomahawk sounding 
rocket 18:63UE. Two concentric deflection plates separated by 0.635 cm 
were held in place by kel-F structures. The inner plate radius , R^, 
was 6.35 cm The plates were in the shape of spherical triangles with 
vertex angles of 90°, 90° and 120° (a central vertex angle of 180° rather 
than 120° would describe a quadrispherical analyzer). Figure 1 is a 
view of the analyzer with the outer plate removed showing the double 
entrance and exit slots as well as typical allowed trajectories. The 
entrance apertures were two 7*37 cm slots on either side of the central 
vertex angle.. These slots were convered by a fine, high transmission 
tungsten mesh grid which was at ground (vehicle) potential. The exit 
apertures were 1.27 cm x 0.635 cm slots in the Kel-F at right angles 
(on the spherical triangle) opposite the entrance slots. Twenty stage 
aluminum dynode electron multipliers were mounted adjacent to the exit 
apertures and were used to detect transmitted electrons. 

Nominal electron trajectories required electrons to enter essentially 
perpendicular to the entrance slots. For a nominal trajectory the upper 
slot was oriented to detect electrons at an angle of 10° from the payload 
axis of symmetry (0° being the upward direction). The nominal lower slot 
viewing angle was 70°, from the spin axis. The minimum arc length traveled 
in a nominal trajectory was tt/ 2 R^. The paths of the nominal electrons 
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from each entrance slot to the designed exit slot actually intersected 
between the plates. During pre-launch calibration there 1 was only time 
for testing the nominal trajectories. After launch subsequent testing 
of an identical set of deflection plates revealed that the exit slot 
collimation provided by the 0.635 cm thick Kel-F was not sufficient 
to exclude some electrons which, not entering normal to the entrance 
slot, actually exited at the nearest exit slot and had a path length 
between the plates less than tt/ 2 R . This degradated the 
response of the analyzer because there was no way of knowing through 
which entrance slot the electron responsible for an anode pulse on the 
electron multiplier had passed. This necessitated assuming that a given 
electron multiplier actually counted electrons which had passed through 
either of two electrostatic analyzers. For example the electron multi- 
plier which nominally accepted electrons viewed at 70° to the spin axis 
had the nominal analyzer with a geometric factor designated as the normal 
geometric factor plus the unintended upward viewing analyzer with a 
geometric factor called the upper slot geometric factor. The total 
analyzer geometric factor was the sum of the response from each slot. 

The upper slot geometric factor had poorer energy resolution than the 
normal geometric factor due to the shorter path length between the 
plates. Even more damaging was the effect of the upper slot 
acceptance upon the angular resolution of the detector. At times 
during the flight of 18:63 UE the pitch angles of electrons arriving 
at the exit slot may have differed by as much as 50°. 



B. Electronics 

The electronics section of the PESFEC provided high voltage for 
the electron multipliers, the positive and negative voltage sweeps for 
the deflection plates and amplified and counted the anode pulses from t 
electron multipliers. The sweep generator and counting section were 
given timing commands from a digital programmer which constructed the 
frame and word intervals from the 10 kHz PCM clock signal. The 0.1088 
second PCM frame consisted of 32 3.4 msec data words . During the first 
3 words of each frame the voltage sweep generator recharged symmetrical 
positive and negative RC networks. The voltages on the capacitors in 
the RC networks decayed with the RC time constant. These voltages were 
used to provide essentially exponentially decaying potentials to the 
deflection plates. Figure 2 shows the deflection plate voltage at the 
start of the accumulation interval for each of the 29 data words. The 
counts in. these -29 data words were used to determine the differential 
electron energy spectrum. The initial three data words did not repre- 
sent an accumulation of counts from the pulse amplifier, hut each con- 
sisted of an identical bit pattern which was used to provide a frame 
sync. 

C. Flat Spectrum - Fixed Plate Voltage Techniques 


Theodcridis and Paolini. define an energy-geometrical factor F 

2 

(cm -sr-keV) as the ratio of the transmitted particle flux (particles/ 

2 

sec) to the directional intensity (particles/sec-cm -sr-ke.V) of ambient 

particle flux. This single parameter relationship assumes that the 

directional intensity is independent of energy over the range of the 

4 

analyzer energy response. Heikkila et al use a similar assumption in 
relating the observed count N.to the differential number spectrum, dj /d) 


_dj_ = N 

dE G * R • E *x *ri 

o UJ 


-2 -1 

cm -sr 


eV 


-1 


sec 


-1 


( 1 ) 


N • (NNF) 





where N is the counts/word, G the geometric factor (era^-sr), R 
the energy resoltuion AE/E q , E q the center energy of the sample 
(in eV), x^ the. sample accumulation time (in sec), n the particle 
detection efficiency and NNF the number normalization factor which is 
very closely related to the reciprocal of the energy-geometrical factor 
F. 

When the differential spectrum is independent of energy (a 

flat spectrum) and the voltages on the analyzer plates are fixed 
during the time x^ , equation (1) is correct and a center energy E q 
can be precisely defined. When the analyzer deflection plate voltages 
are swept rather than stepped from one constant value to a new constant 
value the center energy E q must be defined in terms of the average 
electron energy of the electrons detected during the sample accumulation 
time . Many experiments designed to measure the auroral electron 

energy spectrum have observed not flat spectra but narrow peaks of 
<■ 1 keV width fcee the review by Hones et al~\ Consequently unless the 
energy resolution of the. detector AE/E is less than the change in 
deflection voltage AV/V during the accumulation time x^ there exists 
the very real possibility that most of the accumulated counts may have 
come from electrons with energies closer to the flat spectrum center 
energy of an adjacent data word. Such high resolution detectors are 
generally avoided because the consequent small energy-geometrical factor 
F inhibits the accumulation of significant counts with high time reso- 
lution, In this paper we will describe a differential energy spectrum, 
unfolding technique which uses the accumulated counts from several sequential 
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data words and does not require a flat spectrum assumption. This 

technique is specifically adapted to the case of an electrostatic analyzer 

spectrometer where the plate voltages are exponentially decaying over 

the sample accumulation time interval. Because AE/E 'i 0.4 while 

AV/V h- 0.15 we call the PESPEC analyzer a low resolution analyzer. 

II . Analyzer Respo nse to Isotropic Flux Deflection Plate Voltages 
wit h Constant . 

A. Definition of Variables 

We have adopted the coordinate system used by Theodoridis and 
1 3 

Paolini' and Smith and Day . The angle a , defined in the plane, of 

the trajectory, is the angle of incidence of the electron with respect 

to the normal to the entrance slot. The angle 3 is measured in a 

plane, which is normal to the plane of the entrance slot and tangent to 

the deflection plates at the point of entrance to the plates. Aa is 

the range of the- angle a for trajectories which remain between the 

inner and outer plates (see Figure 3). A3 is the range of the angle 

3 for which the plane of the trajectory is such that the electron can 

pass through the exit slot, (see Figure 4). Because the values of 

Aa and A3 vary with position within the entrance slot we subdivide 

the entrance aperture into many smaller apertures each with some small 

area AA over which Aa and A3 are assumed constant. See Figure 5 

for the definition of the polar coordinates (r » Y.) of the entrance 

i ^ 

slot area AA (note that the indices i and j are decoupled). For 

each entrance slot the geometric factor is the sum of the geometric factors 

of each of these small subdivisions. If the exit aperture were also so 

large that Aa. . or A3., from a given entrance subdivision would not 
i.l ij 

be constant over the area of the exit aperture one would also have to 
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subdivide the exit aperture. The solid angle, , subtended by the 

] 

subdivision AA^ can be determined [see Theodoridis and Paolini , p. 630] 
from 


■Q.’. = 2 A$. . sin(Aa../2) 
13 iJ 13 


( 2 ) 


AA can be determined from the radial, Ar, and angular, Ay , 
spacing between subdivisions and the radial distance, r 


AA. . = (Ar)(Ay)r 
ij o, 


(3) 


For the upper slot where the center of the AB_ range is not 
normal to the AA^ one must use the area projected by the AA__ 
normal to the center value of B . The values of Aa__ depend upon 
the incident electron energy, E , and the plate voltage, V , and there- 
fore the geometric factor for each entrance slot can be written 


i j 

max max 


G(E,V) = l l 

i“l 3=1 


(AA. .)(«..) . 
13 13 


(4) 


The values of the upper geometric factor and the nominal geometric 
factor can be determined by essentially the same method, the only 
difference being that the central angle, (j) , varies with entrance 

°j 

position (r ,y.) for the upper slot. The central angle, <f> , and 

°i J °3 

the polar coordinates (r , 4>) which describe the position of the electron 

between the plates and the plate radii are also illustrated in Figure 3. 
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Neglecting fringing fields the limiting maximum and minimum values 

of a. a and a . , as well as 8 and Af3 can be determined 

’ max min- 

from purely geometrical considerations. 

B. The Trajectory Between the Plates 

For a central electrostatic force one can show (see Paolini and 
2 

Theodoridis ', p. 581) that the trajectory between the plates of the 
electron which enters the plates with initial kinetic energy E and 
angle a is 

r ^ = U(1 4- e cos ({f>— cf> T ) ) (5) 


where is the angle of apsides and the eccentricity, £, is given 

by 


£ “(1 


+ 


2E’ 1/2 
KU ' 


( 6 ) 


where the total electron energy between the plates, E' , is given 
by 


E' 


V (< 0) and V (> 0) 
e is the magnitude of 
K is given by 


E 


r-e(V R -V.R.) 
00 11 

R -R. 


< 0 


(7) 


u 01 

are the outer and inner plate voltages respectively; 
electronic charge. 


-e R R. (V -V.) 
0101 


R -R. 
o x 


K 


( 8 ) 
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| 

The constant U is defined by 


U = 


K 

-12 ■ 2 

■ 2 (E f + K r )r cos a 
o . o . 

X X 


(9) 


The angle of apsides can be determined from the slope of the trajectory 

at r and the derivative of equation (5) with respect to time 
i 


- r ^ r = - U e’sinCcj)-^' )$ 


tan a 


• /A A 1 \ i- 0.1.1. 

SInOM ) - lr5 — 


( 10 ) 


At the entrance (r = r )A = 0 , therefore 

o 

x 


. , _ - tan a 

Sln * " ue 

O . 

X 


(ID 


Also using equation (5) 


-1 

r 

o . 
x 


U(1 + e cos(-cJ>')) 

1-U r 


COS (})' = 


o . 

X 


Ufr 


o . 

X 


( 12 ) 


Equations (11) and (12) uniquely determine <j)' , 


(J)* = tan 


-1 


- tan a 

a-vr'l. 

x -* 


(13) 
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Using the parameters determined in equations (6), (7), (8), (9) 

and (13) one can use equation (5) to determine the trajectory of the 

electron as it passes through the plates for an initial E, r and 

i 

a. 

C. Calculation of Limiting Orbits 

However, an electron may have a trajectory which does not strike 
the plates but is intercepted by the exit collimator, r and a are 
defined to be the exit radius and angle a respectively when <j> = 4> 0 ■ 
After exiting from the plates the electron is assumed to follow a straight 
ray trajectory, can be detex-mined using equation (10) 

tan a = U£r sin(<j> - cj)’) . (14) 

X X Oj 

Defining d(d « R^) to be the distance in the plane of the 
trajectory which the electron must drift to pass the exit collimator 
we can find the radial distance the electron will drift, Ar, by 

' Ar = d tan . (15) 

Therefore the radial position of the electron as it exits the 
collimator, r c , is given by 

r = r + Ar (16) 

c x 

The value of Aa. . is determined from the maximum and minimum 
ij 

values of a within the interval a . < a < a which have 

mm . . max . . 

ij 13 

trajectories which are always between the plates and satisfy < r c < R q 



The computer program which determines Aol. f° r 8i- ven , 

J i 

d) , etc. first determines a and a . . Beginning at ot . 

Y o . ’ max.. min. mlu - ; 

J ^ i , 

test values of a within this interval are used to compute the parameters 

in equations (5) and (16) . Calculating time is minimized by checking 

first for exit clearance and then varying cj> in equation (5) from 

j 

back to 0° in 1° steps. At each step r is calculated to determine 
whether the electron is still between the plates. The first value of 
ot which has an allowed trajectory is defined to be cx^ . ct is increased 
until the trajectory hits the plates or a = . Defining the last 

ij 

allowed trajectory as > Aa_ S^ ven ky 

Aa. . = a 0 - a. . (17) 

ij 2 1 

We emphasize that only minimal computation time (less than 5 
^minutes on a UNIVAC 1108) was required to determine the geometric factor, 
energy response and the allowed orbits calculations used for comparison 
with calibrations. The computations for the nominal slot where the 
energy resolution was not a function of y represented a small 
fraction of this amount. 

D. Comparison with Laboratory Calibrations 

Laboratory measurements of the Aa„ an <3 A3^ were also performed. 

An electron gun which could be varied in a and 3 independently was 

directed at various points (r ,y.) along each entrance slot. Measure- 

°i 

ments x-/ere made at a beam energy of 5 keV because it was sufficient to 

allow the use of phosphorescent screens to determine (r , y.) and 

°i 

insure that the beam diameter was less than 0. 200 cm Uncertainties in 

r Q were 'u 0.10 cm Uncertainty in y_. was 'v 1°. The deflection plate 

voltage, V, was varied rather than electron energy, E, to facilitate 

maintaining constant beam current. Electrons transmitted through the 

plates were collected by a Faraday cup, and the current was measured by an 

electrometer. Another moveable Faraday cup could be positioned to measure 

/ 

the electron gun beam current before it entered the plates. 
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At each position, (r ,Y.) along the slots one could measure the 

°i ^ , 

current transmitted while varying either a , 3 or the deflection plate 
voltage V while keeping the other two parameters fixed. The values of 
a, 3 and V at the half-maximum values of the beam current were used to 


determine Aa. . , A3., and the energy resolution. The 1/4 and 3/4 
ij i.l 

values of the beam current determined the uncertainties. Absolute deter- 


minations of a and 3 were difficult because they required a very 


accurate positioning and measurement of orientation of a very irregularly 

shaped detector in the cramped quarters of the vacuum chamber. 

Because a comparison with theoretical determinations of Aa.. 

il 

required the absolute value of a we followed the practice of Paolini 

2 

and Theodoridis of translating the data points in a-space to give the 
best fit between the theoretical and experimental limiting values of 


a and V. The maximum shift needed was 3.3°. The points were not 


shifted in V-space. 


Figure 6 depicts the limiting values in (a,V) space for the nominal 


slot with r = 6.67 cm and Y. = 45°. The shift in absolute a was 
o . j 

l J 

-2.4°. The computed limiting values are for r = 6.76 cm. This was 

i 

the closest theoretical value of r to the experimental value vriiich 

o . 

l 

was computed. The experimental value of A3. . for this slot with 


y. = 45° was 
J 

A3. . = 7.2° ± 2.7° 
ij 


compared to the theoretical value of 5.2°. 


\ 
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Figure 7 shows the experimental and computed limiting values in 

(a,V) space for the upper entrance slot with y_. = 45°. The experimental 

value of r was 6,59 cm and the nearest computed value was 6.63 cm. 
i 

Note that both the experimental and computed limiting values show the 
effect of collimation upon the negative a - small V boundary. Also the 
Aa-AE resolution is much poorer for the upper entrance slot because the 
central angle <|> is less than 90°. For this slot and position the. 
experimental value of Af3„ was a l so consistent with the theoretical 
value. 

Figure 8 shows the geometric factor from each slot and the total 
geometric factor as a function of E/V where V is the plate voltage for 
symmetrical positive and negative voltages on the plates. The actual 
plate voltages on the 18:63 PESPEC deviated somewhat from symmetry due 
to the difference in capacitance of the inner and outer plate plates 
themselves, but the computer results demonstrated that for the small 
departures from symmetry on the PESPEC the effects were negligible. 

The angular resolution as well as the energy resolution of the 
18:63 PESPEC was impaired by the electrons which were able to enter the 
upper slot and exit at the wrong exit operture. The angular response 
of the nominal entrance slot was n, 6° x 6°. The acceptance direction 
in a plane perpendicular to the spin axis (look azimuth) was 123.5° 
from the payload x-axis reference. Because this slot accepted particles 
70° from the spin axis the look elevation was .20°. The upper slot 
acceptance direction was determined by analyzing the slot response in 


terms of the response of five separate detectors with 


values of 19.5°, 



28.5°, 37.5°, 46.5° and 55.5°. By weighting the acceptance azimuth and 
elevation of each of the five detectors with its appropriate energy-angle 
factor an effective acceptance azimuth and elevation could be determined. 
Table J lists the values of the azimuth and elevation for each slot and 
the approximate angular resolution. 

Because 18:63 UE went into a near flat spin attitude during despin the 
difference in the pitch angles of the electrons transmitted through each 
slot varied from 0° to ^ 55° twice per roll. Analysis of pitch angle 
information must be restricted to those portions of each roll where there 
is some equality between the pitch angles of the electrons transmitted 
through the separate slots. 

III. Analyt ic Energy Spect ru m Unfolding Te c hnique 
A. Mathematical Derivation 

We want to determine the differential energy spectrum, dj(E)/dE 
-2 “1 1 

[electrons-sec -cm -sr -KeV ] from a system of equations of the form 


N. 

l 



) n(E) G (E/V (t) ) 
J 0 


dj(E) 

dE 


dE 


(18) 


The 29 data words measured during the voltage sweep on the PESPEC 
give the counts 1L. is the sample time (0.0032sec). rq (E) is the 

electron multiplier detection efficiency as a function of energy E. 

G is the geometric factor shown in figure 8. We assume that the are 

corrected for dead time counting losses. 

Because W(E) is a weak function of energy and the function G(E/V(t)) 
effectively samples only a narrow range of energies one can compute an 
average efficiency for data word i, r| , from 
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T <» 

t 0) 


n. = 


_ •'o 


dt 


n(E) G(E/V(t))dE 


1 T 

r w 


dt 


G(E/V(t))dE 


(19) 


- n(E_ ) 


where E is the center energy of data word i . 
o 

l 

Using n. equation (18) can be written 


N. 

l 


= n. 


0) 


dt 


G (E/V(t) ) - dE 


- (20) 


We assume that we can express dj(E)/dE as a polynomial of order 

(j -1) in E, 

J inax ’ 


dj (E) 
dE 


max 


= l C,E j_1 . 
j=l 3 


(21) 


Define I. by 
J 


r 


I. = 

3 


G (E/V(t) ) E j_1 dE . 


( 22 ) 


Equation (20) can now be written as 


n. = n. 

i i 


T j 
t oj max 

dt }' C.I. 

0 A j j 


(23) 


Because G is a function of the ratio E/V(t) numerical integration 


of equation (22) gives 
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I. = F (t) (24) 

where the values of the constants F. for j = are given 

in table II. 

Interchanging the order of integration and summation in equation 
(23) gives 


N. 

l 


i J 


max 


T 


3 = 1 


C.F. 
J 3 


V (t)*^ dt . 


(25) 


Over the accumulation time interval the PESPEC deflection plate 

voltage can be very well represented by an exponential decay with time 

constant T. from the value V 

1 o. 


V(t) = V e 
o. 

l 


-t/Ti 


(26) 


Using equation (26) the integral in equation (25) can be evaluated 
tp give 


max I -| 

N. = O.T. I C.F. — 

1 1 1 j=l 3 J L 


-jT /T." 
J 0) 1 


o . 

1 


(27) 


Defining 


F. 

H. . 5 I 

j L 


1-e 


-jT / T. 
J lx)' 1 


]• 


(28) 


a count rate, R. , 
’ l ’ 


N. 


R. E — , 

1 U. T. V ' 

i i o . 

l 


( 29 ) 
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and factoring out V from equation (26) gives 

i 

j 


max . , 

R. = l C. H. . V J 
1 3=1 3 °i 


(30) 


Equation (30) describes a system of simultaneous equations which can 
be solved for the values of (h allowing one to piecewise determine 
dj (E)/dE. 

B. Method of Application to Low Resolution Detector 
Exact solutions to the system of equations (30) for higher 
order polynomials (j > 3) may display erratic behavior between the 
fitted points. Closer examination of equation (30) reveals that defining 






(31) 


one can express R. as a polynomial and compute the using a least 

•t J 

squares fit. The approximation in equation (31) would be exact if the 
complete voltage sweep could be fitted by a single decay time (see equation 
(28)). 

Table III lists the 18:63 UE PESPEC values of V , T., r|, and 

o . i i 

l 

H. . for j =5. The values for V , T. and r) . were determined 
i j max o . i i 

i 

using prelaunch calibrations. 

Table III shows that for five point fits the values of Tg are 
constant to ^3%, except near the beginning and end of the sweep where 
over five data words the vary by 15%. 

A weighted least squares fit applied piecewise to the system of 


polynomials 
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max 


R. = ) B. V 3 

i J o. 

1=1 J 1 


3-1 

j 


(32) 


gives the coeff icieiits C^. of equation (21) from equation (31) 


C. = B'JW. . . 
J 3 iJ . 


(33) 


Of course these coefficients are valid only for some energy interval 
near the center energy of the middle word of the piecewise least squares 
fit. 

For a flat spectrum, the average or center energy measured during word 

i, E , can be determined by first numerically computing equation (22) 
i 

for j = 1 . 


I = I G(E/V(t))dE = F V ( t ) . 

'o 


(34) 


then E is defined by 
°i 


E 

o . 

1 1 
7 F V(t) = G(E/V(t))dE 

2 1 0 


(35) 


Equation (35) can be numerically solved for E 

i 


E = 10.6 V (36) 

o . 

l 

where V is the average value of V(t) during the time T . Evaluating 
V we obtain 



10.6 V T. 


o . 1 
i 


T 


to 


fl 



(37) 
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Table IV lists the E values of the PESPEC for each of the 29 data 

°i 

words per frame. 

A correction to the E due to a non-zero slope, S , in the 

i 

energy spectrum can be derived. For a square box geometric factor with 
energy resolution AE/E = 26 and a first order energy dependence for 
it can be shown that the average energy of the electrons detected 
during word i, E_^ , is 1 

5 2 + 0(6 4 ) . (38) 

7 

2 

The best value of 6 for the. actual PESPEC geometric factor was 

numerically computed for realistic values of S , E and dj (E )/uE. 

2 0± 0± 

The best value of 6 depends upon the sign of the slope. This is 

reasonable in view of the skewed geometric factor. For a positive slope 

we find 

6 2 = 0.12 

and for a negative slope we obtain 

S 2 = 0.07 . 

The proper weights to use for the piecewise fitting of the system 

of equations (32) can be determined by computing the fraction of the 

counts measured in the middle word of the fit, i Q , which are actually 

due to electrons with energies nearer the center energies, E , of 

i 

adjacent data words. Figure 9 shows for a flat energy spectrum the 
fraction of the counts accumulated in the middle word which are due to 
electrons with energies closer to the center energies of adjacent words. 


<W 

i 


= 1 + 


Je (E o.>> 

i 
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Note that these fractions which are used as the weights for the fitting 

procedure peak at the middle word i . For a flat spectrum only 'u 38% 

of the counts are due to electrons in the i^ word energy band; however 

when a five point fit is used over 85% of the counts are due to electrons 

with energies in the range of the fitted points. Also shown in figure 9 

are the weights which would be used to unfold the spectrum coefficients 

if only the high resolution normal geometric factor were used. 

The algorithm for the piecewise unfolding of the energy spectrum 

is complete. If for example we choose to fit the counts from 5 words 

for up to second order energy dependence in dj(E)/dE we begin with word 

3 and fit equation (32) from i = 1 to i = 5 with j =3. The 

max 

values of C. are determined from 
3 


C. - . 

3 3 3j 


Using these C. and E = E from equation (37) d j (E )/dE from 
3 ° 3 °3 

equation (21) is computed. The slope at E = E can also be computed 

°3 

using the . Inserting these parameters into equation (38) the corrected 

average energy, E^ , of the electrons counted during word 3 is computed. 

Inserting this corrected average energy into equation (21) gives the value 

of the differential flux at E^. This process is repeated for word 4 

except that the fit begins with counts from word 2 rather than word 1. 

This procedure is then repeated through word 27. 

It is interesting to evaluate this procedure for j ^ = 1 in the 

limit T /t. « 1 which would correspond to constant deflection plate 
Cl) i 

voltages and a histogram type electron differential energy spectrum. 
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therefore 


-T / T. 

[1 - e U 1 ] ->- T /x. 

(O X 


E 

o . 

1 


+ 10.6 V 


o . 

i 


H... +F. T /t. 

il 1 (o' l 


i — -*■ C.F.x /t. 

n.x.v l l to i 


110 , 


and 


N. 


1 F, »V • t *ri. 

1 o . U) 1 
i 


as one would expect from equation (1) . 


IV . Examples of Spectrum Unfolding Technique 

Parameterization of the auroral electron, energy spectrum allows 
one to numerically compute the in equation (18) for various shapes 

of the energy spectrum. By applying the unfolding procedure to the 
computed one can determine how well the unfolded spectrum matches 

the input spectrum. This allows one to determine the optimum j 
and number of values of over which equation (32) should be fitted. 

One can also compare the unfolded spectrum to the spectrum determined 
using equation (1) by unfolding with j = 1. 

2 

The auroral electron differential energy spectrum [electrons/cm - 
sec^-sr-keV] can be represented by an equation of the form 
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f|(E) = J E " L + J E " H + 0.2X10 10 • ■ E 

L H e 


EXP (-(E + E - 2/E£„) /T ) 
D D e 


(39) 


J =0 for E > E 

°l ' C 

J = 0 for E < E 

°H C 


J n T . J and n., are the usual power lav; parameters for energies 
°L L ' °H 11 

below and above E„ respectively, n and T are the directional density 

C e e 

and temperature of a Maxwellian electron gas drifting with respect to the 
detector. The parameter E^ is the equivalent kinetic energy of an electron 
moving at the drift velocity. The peak in the energy spectrum is due to 
the drifting Maxwellian electrons. The width of the peak is proportional 
to the parameter T e . 

Using representative values of the parameters in equation (39) 
we have found that three is the optimal value of the parameter j t0 

use in piecewise unfolding PESPEC spectra. This would correspond to a 
quadratic energy dependence of the energy spectrum over the energy reso- 
lution of the analyzer. Similarly the optimal number of Ik over which 
the piecewise least squares fit is to be performed was found to be five. 

This is consistent with the range of significant weights in figure 9. 

Figure 10 shows the values of Ik which would be measured for two 
energy spectra which represent typical values of the parameter T^ . 

Table V lists the values of the parameters in equation (34) for these hot and 


cold energy spectra. 
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Figures 11 and 12 show the actual values of the differential 

flux, the five point, j = 3 unfolded values and the single point, 

''max = ^ s P ectrum which would be obtained using equation (1) . Figures 

13 and 14 show the fractional error between the unfolded spectra and 

the actual spectrum. Hotter spectra can be unfolded more accurately 

than those with 0.10 keV. The average percentage error for the 

1 =3 unfolded hot spectrum was ^ 3% while it was 'V 5% for a 

max 

= 0.090 KeV spectrum. In the peaked region (data words 4 through 15) 
the . T = 1.26 KeV spectrum has a percentage error of ^3%, but the 
percentage error for the more difficult to unfold = 0. 330 keV spec- 
trum is n. 7%. We note that the j ^ = 3 unfolding has preserved 
the width of the peak and the value of dj/dE at the peak much better 
than the j =1 technique (see figures 13 and 14) . 
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Table I. Azimuth and Elevation of each Detector Entrance Slot 


and Approximate Angular 

Resolution 


Slot 

Azimuth 

Elevation 

Angular Resolution 

Nominal Slot 

123.5° 

O 

O 

CM 

-6° x 6° 

Upper Slot 

180.8° 

42.4° 

~7° x 35° 
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Table II. Values of Constants Determined from Numerical 
Integration of Equation (22). 


j 

1 

2 

3 

4 

5 



0.223 

0.226 x 10 1 
0.360 x 10 2 
0.584 x 10 3 
0.119 x 10 5 




) 
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Table III. 

18:63 

UE PESPEC 

Voltage 

Sweep 

and Analyzer 

Parameters 

i 

V [kV] 
o . 
i 

T_^[sec]ru 

H il 

H i2 

H i3 

H i4 

H i5 

1 

3.933 

.0164 

.100 

.0395 

.430 

5.31 

79.1 

1480. 

2 

3.197 

.0166 

.100 

.0392 

.426 

5.28 

78.6 

1470. 

3 

2.604 

.0171 

.100 

.0381 

.415 

5.16 

76.9 

1450. 

4 

2.135 

.0195 

.100 

.0338 

.372 

4.67 

70.3 

1330. 

5 

1.793 

.0212 

.100 

.0313 

.347 

4.37 

66.2 

1260. 

6 

1.527 

.0214 

.100 

.0310 

.344 

4.34 

65.7 

1250. 

7 

1.303 

.0212 

.100 

.0312 

.346 

4.37 

66.1 

1260. 

8 

1.110 

.0215 

.100 

.0308 

.342 

4.32 

65.5 

1250. 

9 

.948 

.0217 

.108 

.0306 

.340 

4.30 

65.2 

1240. 

10 

.810 

.0221 

.124 

.0301 

.335 

4.23 

64.3 

1230. 

11 

.694 

.0220 

.138 

.0302 

.336 

4.24 

64.4 

1230. 

12 

.595 . 

.0225 

.152 

.0296 

.330 

4.17 

63.4 

1210. 

13 

.511 

.0226 

.165 

.0294 

.328 

4.15 

63.1 

1210. 

14 

.440 

.0227 

.177 

.0294 

.327 

4.14 

63.0 

1200. 

15 

.379 

.0226 

.188 

.0295 

.329 

4.16 

63.2 

1210. 

| 

16 

.326 

.0229 

.198 

.0290 

.324 

4.10 

62.4 

1200. 

17 

.281 

.0232 

.207 

.0287 

.321 

4.06 

61.9 

1190. 

18 

.243 

.0236 

.215 

.0282 

.315 

4.01 

61.0 

1170. 

19 

.210 

.0240 

.223 

.0278 

.311 

3.96 

60.3 

1160. 

20 

.182 

.0240 

.230 

.0278 

.311 

3.95 

60.3 

1160. 

21 

.158 

.0244 

.238 

.0274 

.307 

3.91 

59.7 

1150. 

22 

.138 

.0249 

.245 

.0269 

.302 

3.84 

58.8 

1130. 

23 

.120 

.0252 

.252 

.0266 

.299 

3.81 

58.2 

1120. 

24 

.105 

.0257 

.258 

.0261 

.293 

3.74 

57.3 

1100. 

25 

.092 

.0262 

.265 

.0256 

.288 

3.68 

56.4 

1090. 

26 

.081 

.0268 

.272 

.0251 

.282 

3.61 

55.4 

1070. 

27 

.071 

.0303 

.277 

.0224 

.253 

3.26 

50.3 

977. 

28 

. .064 

.0279 

.283 

.0242 

.273 

3.50 

53.8 

1040. 

29 

.056 

.0347 

.288 

.0196 

.224 

2.90 

45.0 

879. 
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Table IV. Spectrometer 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


Center Energies for 

E [keV] 
o . 

l 

37.8 

30.8 
25.1 

20.8 
17.6 
15.0 
12.8 
10.9 

9.33 
7.99 
6.84 
5.87 
5.05 

4.34 
3.74 
3.22 
2.78 
2.40 
2.08 
1.80 
1.57 
1.36 
1.19 
1.04 

.916 
.806 
.715 
.636 
.569 


Flat Energy Spectrum. 


> 



Table V. Electron Energy Spectrum Parameters for Spectra 
Used to Illustrate Spectrum Unfolding. 


Spectrum 

J, 

°L 

7 

n 

L 

J 

°H 

7 

n 

H 

Hot 

5.24x10 

0.593 

3.63x10 

0.504 

Cold 

2. 14xl0 7 

0.583 

7. 11x10 8 

1.768 


r -3 -1, 

n [cm -sr J 
e 

-3 

T [keV] 
e 

E D [keV] 

E c [keV] 

1.36x10 

1.260 

11.0 

28.42 

-L 

1.18x10 

0.330 

11.4 

16.75 
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FIGURE LEGENDS 


Figure 1. 


Figure 2. 


Figure 3. 


Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


View of analyzer with outer plate removed to show double 
entrance and exit slots. Nominal trajectories as well as 
those from the upper slot are indicated. The analyzer is 
resting on the electron multiplier. 

Magnitude of the voltage on deflection plates at start of 
the counting interval for each of the 29 data words. Note 
the approximately exponential decay character of the voltage 
sweep . 

View in the plane of trajectory of electron as it passes 
between the plates illustrating the coordinates a, r, $ 
and the analyzer geometrical parameters. 

View of the analyzer in the plane of the Kel-F baseplate 
containing exit slots. The parameters 3 and A3 are shown. 
View of an entrance slot illustrating the subdivisions and 
the coordinates of the subdivisions. 

Experimental and calculated limits of allowed orbits in 
(a,V) space for normal entrance slot and 5 keV electron 
energy. 

Experimental and calculated limits of allowed orbits in 
(a,V) space for upper entrance slot and 5 keV electron 
energy. Note change of abcissa scale from fig. 6. 



31 


Figure 8. 


Figure 9. 


Figure 10. 


Figure 11. 


Figure 12. 


Figure 13. 


Figure 14. 


Geometric factor for each slot and total geometric factor 
versus electron energy when deflection plate voltages are 
± 10 volts. 

Weights used in piecewise fitting systems of equations (32) 

about data word i . 

o 

Counts per data word, N^, resulting from the two differential 
energy spectra described in Table V. When the data word 
numbers are in reverse order the abcissa is very nearly a 
logarithmic energy scale. 

Actual differential energy spectrum and the j = 1 and 

^max = ^ un f°l<ied spectra in the energy range of the peak. 

The T g = 1.26 keV spectrum of Table V is the actual spectrum. 

Actual differential energy spectrum and the 3 max = 1 and 

j max = 1 unfolded spectra in the energy range of the peak. 

The T & = 0.330 keV spectrum of Table V is the actual spectrum. 

Note that the j ^ = 1 unfolded spectrum gives a much more 

accurate description of the width and height of the peak. 

Fractional error between unfolded spectra and the actual 

spectrum over complete detector energy range for T =1.26 keV 

e 

^ j 

spectrum. 

Fractional error between unfolded spectra and the actual 
spectrum over complete detector energy range for T^ = 0.330 keV 
spectrum. 






























* 6.67cm EXPERIMENTAL 
‘i 6.67cm CALCULATED 



NORMAL ENTRANCE SLOT 
ELECTRON ENERGY = 5.0 keV 

VOLTS 


















COUNTS N i 






o 



II 






FRACTIONAL ERROR 




FRACTIONAL ERROR 






